Effective Stiffness: Generalizing Effective Resistance Sampling to Finite Element Matrices
نویسندگان
چکیده
We define the notion of effective stiffness and show that it can used to build sparsifiers, algorithms that sparsify linear systems arising from finite-element discretizations of PDEs. In particular, we show that sampling O(n log n) elements according to probabilities derived from effective stiffnesses yields an high quality preconditioner that can be used to solve the linear system in a small number of iterations. Effective stiffness generalizes the notion of effective resistance, a key ingredient of recent progress in developing nearly linear symmetric diagonally dominant (SDD) linear solvers. Solving finite elements problems is of considerably more interest than the solution of SDD linear systems, since the finite element method is frequently used to numerically solve PDEs arising in scientific and engineering applications. Unlike SDD systems, which are relatively easy to precondition, there has been limited success in designing fast solvers for finite element systems, and previous algorithms usually target discretization of limited class of PDEs like scalar elliptic or 2D trusses. Our sparsifier is general; it applies to a wide range of finite-element discretizations. A sparsifier does not constitute a complete linear solver. To construct a solver, one needs additional components (e.g., an efficient elimination or multilevel scheme for the sparsified system). Still, sparsifiers have been a critical tools in efficient SDD solvers, and we believe that our sparsifier will become a key ingredient in future finite-element solvers.
منابع مشابه
Nonlinear inelastic static analysis of plane frames with numerically generated tangent stiffness matrices
For the nonlinear analysis of structures using the well known Newton-Raphson Method, the tangent stiffness matrices of the elements must be constructed in each iteration. Due to the high expense required to find the exact tangent stiffness matrices, researchers have developed novel innovations into the Newton-Raphson method to reduce the cost and time required by the analysis. In this paper, a ...
متن کاملA GUIDED TABU SEARCH FOR PROFILE OPTIMIZATION OF FINITE ELEMENT MODELS
In this paper a Guided Tabu Search (GTS) is utilized for optimal nodal ordering of finite element models (FEMs) leading to small profile for the stiffness matrices of the models. The search strategy is accelerated and a graph-theoretical approach is used as guidance. The method is evaluated by minimization of graph matrices pattern equivalent to stiffness matrices of finite element models. Comp...
متن کاملA homogenization sampling procedure for calculating trabecular bone effective stiffness and tissue level stress.
A homogenization sampling procedure is introduced which allows computation of effective trabecular bone stiffness and individual trabecula level stress based on precise models of trabecular bone architecture. Three-dimensional digitized images of 53 trabecular bone specimens with a resolution of 50 microns per voxel were directly converted into three-dimensional finite element meshes by making ...
متن کاملBuckling Behaviors of Symmetric and Antisymmetric Functionally Graded Beams
The present study investigates buckling characteristics of both nonlinear symmetric power and sigmoid functionally graded (FG) beams. The volume fractions of metal and ceramic are assumed to be distributed through a beam thickness by the sigmoid-law distribution (S-FGM), and the symmetric power function (SP-FGM). These functions have smooth variation of properties across the boundary rather tha...
متن کاملApplication of Decoupled Scaled Boundary Finite Element Method to Solve Eigenvalue Helmholtz Problems (Research Note)
A novel element with arbitrary domain shape by using decoupled scaled boundary finite element (DSBFEM) is proposed for eigenvalue analysis of 2D vibrating rods with different boundary conditions. Within the proposed element scheme, the mode shapes of vibrating rods with variable boundary conditions are modelled and results are plotted. All possible conditions for the rods ends are incorporated ...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
عنوان ژورنال:
- CoRR
دوره abs/1110.4437 شماره
صفحات -
تاریخ انتشار 2011